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ABSTRACT 

The flat spectrum radio quasar (FSRQ) PKS 2123—463 was associated in the first Fermi- 
Large Area Telescope (LAT) source catalogue with the y-ray source 1FGL J2126. 1—4603, 
but when considering the full first two years of Fermi observations, no y-ray source at a 
position consistent with this FSRQ was detected, and thus PKS 2123—463 was not reported 
in the second Fermi - LAT source catalogue. On 2011 December 14 a y-ray source positionally 
consistent with PKS 2123—463 was detected in flaring activity by Fermi- LAT. This activity 
triggered radio-to-X-ray observations by the Swift, Gamma-ray Optical/Near-Infrared Detector 
(GROND), Australia Telescope Compact Array (ATCA), Ceduna and Seven Dishes Karoo 
Array Telescope (KAT-7) observatories. Results of the localization of the y-ray source over 41 
months of Fermi - LAT operation are reported here in conjunction with the results of the analysis 
of radio, optical, ultraviolet (UV) and X-ray data collected soon after the y-ray flare. The strict 
spatial association with the lower energy counterpart together with a simultaneous increase 
of the activity in optical, UV, X-ray and y-ray bands led to a firm identification of the y-ray 
source with PKS 2123—463. A new photometric redshift has been estimated as z — 1.46 ± 
0.05 using GROND and Swift Ultraviolet/Optical Telescope (UVOT) observations, in rough 
agreement with the disputed spectroscopic redshift of z — 1.67. We fit the broad-band spectral 
energy distribution with a synchrotron/external Compton model. We find that a thermal disc 
component is necessary to explain the optical/UV emission detected by SwiftfUVOT. This disc 
has a luminosity of ~1 .8 x 10 46 erg s _1 , and a fit to the disc emission assuming a Schwarzschild 
(i.e. non-rotating) black hole gives a mass of ~2 x 10 9 Mq . This is the first black hole mass 
estimate for this source. 

Key words: galaxies; active - galaxies: nuclei - quasars: general - quasars: individual: 

PKS 2123—463 - gamma rays: general. 


1 INTRODUCTION 

Blazars constitute the most extreme subclass of active galactic nu- 
*E-mail: filippo.dammando@flsica.unipg.it clei (AGN), characterized by the emission of strong non-thermal 
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radiation across the entire electromagnetic spectrum and in partic- 
ular intense y-ray emission above 100 Me V. The typical observa- 
tional properties of blazars include irregular, rapid, high-amplitude 
variability, radio-core dominance, apparent superluminal motion, 
a flat radio spectrum, and high and variable polarization at radio 
and optical frequencies. These features are interpreted as resulting 
from the emission of high-energy particles accelerated within a rel- 
ativistic jet closely aligned with our line of sight and launched in 
the vicinity of the supermassive black hole harboured by the active 
galaxy (Blandford & Rees 1978; Urry & Padovani 1995). 

Since the advent of the Energetic Gamma-Ray Experiment 
Telescope (EGRET) on the Compton Gamma-Ray Observatory 
(CGRO), blazars were known to dominate the extragalactic high- 
energy sky. However, EGRET did not pinpoint the location of many 
sources with sufficient precision to enable astronomers to associate 
them with known objects, leaving the legacy of a large fraction 
of unidentified sources in y rays. The point spread function (PSF) 
and sensitivity of the Large Area Telescope (LAT) on-board Fermi 
provide an unprecedented angular resolution at high energies for 
localizing a large number of newly found y -ray-emitting sources. 
Correlated variability observed at different frequencies can give im- 
portant information for the identification of a y-ray source with its 
low-energy counterpart. 

On 2011 December 14 a y-ray flare from a source positionally 
consistent with PKS 2123—463 was detected by Fermi - LAT 
(Orienti & D’ Ammando 2011), triggering Gamma-ray 
Optical/Near-Infrared Detector (GROND) and Swift follow- 
up observations (D’ Ammando, Orienti & Mountford 2011; Rau 
et al. 2011) that confirmed contemporaneous activity in the 
optical/ultraviolet (UV) as well as marginally in X-rays. 

PKS 2123—463 is a bright radio quasar with a luminosity at 
1.4GHz Li. 4 ghz = (1.5 ± 0.2) x 10 28 \VHz~ 1 (assuming the red- 
shift estimated in this paper, z = 1.46 ± 0.05 ;' see Section 6). 
On the basis of its spectral index a r ~ 0.4 ()>(v) oc v - “ r ) be- 
tween 408 MHz and 4.8 GHz, it was included in the Combined 
Radio All Sky Targeted Eight GHz Survey (CRATES) catalogue 
of flat spectrum objects (Healey et al. 2007). Australia Telescope 
Compact Array (ATCA) observations performed almost simultane- 
ously at 4.8, 8.6 and 20 GHz during the Australia Telescope 20 GHz 
(AT20G) survey indicated a flattening of the radio spectrum at high 
frequencies to a spectral index of a r ~ 0.2. Polarized emission 
has been detected only at 4.8 GHz where the polarization is about 
4 per cent of the total intensity flux density (Massardi et al. 2008). A 
Chandra observation of PKS 2123—463 in 2004 March has shown 
the presence of a jet-like extended structure in X-rays (Marshall 
et al. 2011). The redshift of PKS 2123—463 was reported to be 
Z = 1.67 (Savage & Wright 1981) based on an objective-prism 
spectrum, subsequently questioned by Jackson et al. (2002) be- 
cause two possible redshifts (0.48 and 1.67) were given from ob- 
servation of two lines in the spectrum, and the motivation for the 
exclusion of the smaller redshift was not provided. This object 
was a member of the pr e-Fermi launch Roma-BZCAT (Massaro 
et al. 2009) catalogue listing candidate y-ray blazars detectable by 
Fermi - LAT but not in the Candidate Gammay-Ray Blazar Survey 
(CGRaBS) catalogue (Healey et al. 2008). 

The flat spectrum radio quasar (FSRQ) PKS 2123—463 (RA: 
21 h 26 m 3K7042, Dec.: -46 d 05 m 47 s .892, J2000; Fey et al. 2006) was 
associated in the First Fermi - LAT source catalogue (1FGL, 2008 


1 The corresponding luminosity distance forr= 1.46 ± 0.05 is^L = 10.6 ± 

0.5 Gpc and 1 arcsec corresponds to a projected distance of 8.53 ± 0.02 kpc. 


August 4-2009 July 4; Abdo et al. 2010a) with the y-ray source 
1FGL J2126. 1—4603, while no association with a y-ray source 
was present in the second Fermi - LAT source catalogue (2FGL, 
2008 August 4-2010 August 4; Nolan et al. 2012), although a 
y-ray source, 2FGL J2125.0— 4632, with a radius of the 95 per cent 
source location confidence region of 0?17, at (L52 from the radio 
position of PKS 2123—463 was reported. Taking into consideration 
the high variability of blazars, the flux of PKS 2123—463 could have 
decreased in the second year of Fermi - LAT operation. In addition, 
in the construction of the 2FGL catalogue, PKS 2123—463 was split 
into more than one candidate source seed (see section 4.2 of Nolan 
et al. 2012, for details). These two factors may have led to the lack 
of its association with a y-ray source. 

In this paper, we present the localization over 41 months of 
Fermi-hAI data of a y-ray source associated with the FSRQ PKS 
2123—463. The correlated variability observed in optical, UV, X-ray 
and y rays confirms the identification. In addition, a new estimation 
of the redshift of the source by means of the fit of simultaneous 
GROND and Swift Ultraviolet/Optical Telescope (UVOT) data col- 
lected soon after the y-ray flare is presented. The paper is organized 
as follows. In Section 2, we report the LAT data analyses and results. 
In Sections 3 and 4, we report the results of the Swift and GROND 
data analysis, respectively. Radio data collected by the ATCA, Ce- 
duna and Seven Dishes Karoo Array Telescope (KAT-7) telescopes 
are presented in Section 5. Section 6 presents an estimation of the 
redshift of the source. In Section 7, we discuss the modelling of 
the overall spectral energy distribution (SED) and draw our conclu- 
sions. Throughout the paper a A cold dark matter cosmology with 
Ho = 71 km s _1 Mpc -1 , £2 a = 0.73 and £2 m = 0.27 is adopted. 

2 Fermi-L AT DATA: SELECTION AND 
ANALYSIS 

The Fermi-CAtY is a y-ray telescope operating from 20MeV 
to >300 GeV. It has a large peak effective area (~8000cm 2 for 
1 GeV photons), an energy resolution of typically ~10 per cent and 
a field of view of about 2.4 sr with an angular resolution (68 per cent 
containment angle) better than 1° for energies above 1 GeV. Further 
details about the Fermi - LAT are given by Atwood et al. (2009). 

The Fermi-LAtT data reported in this paper were collected over 
the first 41 months of Fermi operation, from 2008 August 4 (MJD 
54682) to 2012 January 4 (MJD 55930). During this time the LAT 
instrument operated almost entirely in survey mode. The spectral 
analysis was performed with the science tools software package 
version v9r23p 1 . The Fermi-FPCF data were extracted from a circular 
region of interest (Rol) with a 1 5° radius centred at the radio location 
of PKS 2123—463. Only events belonging to the ‘source’ class were 
used. The time intervals when the rocking angle of the LAT was 
greater than 52° were rejected. In addition, a cut on the zenith 
angle (<100°) was also applied to reduce contamination from the 
Earth limb y rays, which are produced by cosmic rays interacting 
with the upper atmosphere. The spectral analyses (from which we 
derived spectral fits and photon fluxes) were performed with the 
post-launch instrument response functions (IRFs) P7S0URCE_V6 
using an unbinned maximum likelihood method implemented in 
the Science tool gtlike. 

The background model used to extract the y-ray signal includes 
a Galactic diffuse emission component and an isotropic compo- 
nent. The model that we adopted for the Galactic component, 
the same as used for the 2FGL catalogue, is given by the file 
gal_2yearp7v6_v0.fits, and the isotropic component, which is the 
sum of the extragalactic diffuse emission and the residual charged 
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particle background, is parametrized by the file iso_p7v6source.txt. 2 * * 
The normalizations of both components in the background model 
were allowed to vary freely during the spectral point fitting. 

We examine the significance of the y-ray signal from the sources 
using the test statistic (TS) based on the likelihood ratio test. The 
TS, defined as TS = 2Alog(likelihood) between models with and 
without the source, is a measure of the probability of having a 
y-ray source at the localization specified, which compares models 
whose parameters have been adjusted to maximize the likelihood 
of the data given the model (Mattox et al. 1996). The source model 
used in gtlike includes all the point sources from the 2FGL that 
fall within 20° of PKS 2123—463. The spectra of those sources 
were parametrized by power-law (PL) functions, except for 2FGL 
J2056.2— 4715 for which we used a log-parabola (LP), and 2FGL 
J2124.6— 3357 and 2FGL J2241.7— 5236 for which we used an ex- 
ponentially cut-off PL in their spectral modelling as in the 2FGL 
catalogue. We removed from the model the sources having TS < 25 
and/or fluxes (0.1-100GeV) below 1.0 x 10~ 8 photons cm~ 2 s _1 
over 41 months and repeated the fit. We tested whether two distinct 
y-ray sources (one at the radio position of PKS 2123—463 and one 
at the y-ray position of 2FGL J2125. 0—4632) are detected simul- 
taneously by Fermi with TS > 25 over 41 months of observations. 
In this scenario, the fit yields a TS of 553 for PKS 2123—463 and 
a TS of 14 for 2FGL J2125.0-4632. We thus conclude that PKS 
2123—463 and 2FGL J2125.0— 4632 are the same source and we 
maintain only PKS 2123—463 in the model. Thus a final fitting 
procedure has been performed with all sources within 10° of PKS 
2123—463 included with the normalization factors and the photon 
indices left as free parameters. For the sources located between 10° 
and 15° we kept the normalization and the photon index fixed to the 
values obtained in the previous fitting procedure. The Rol model 
includes also sources falling between 15° and 20° from the target 
source, which can contribute to the total counts observed in the 
Rol due to the energy-dependent size of the PSF of the instrument. 
For these additional sources, normalizations and indices were fixed 
to the values of the 2FGL catalogue during all steps of the fitting 
procedure. 

The y-ray point source localization using the gtfindsrc tool 
over the photons extracted during the period 2008 August 4-2012 
January 4 results in RA = 321?609, Dec. = — 46?076 (J2000), with 
a 95 per cent error circle radius of 0?047, at an angular separation of 
0?024 from the radio position of PKS 2123—463 (RA = 321?628, 
Dec. = — 461097, J2000). This implies a strict spatial association 
with PKS 2123—463. The fit with a PL model, dN/dE oc (E/E 0 )~ r , 
to the data integrated over 41 months of Fermi operation in the 
0. 1-100 GeV energy range results in a TS = 557, with an integrated 
average flux of (4.76 ± 0.35) x 10~ 8 photons cm~ 2 s _1 and a photon 
index T = 2.45 ± 0.05 (with E 0 fixed to 670 Me V). In order to test 
for curvature in the y-ray spectrum of PKS 2123—463 we used as 
an alternative spectral model the LP, dN/dE oc (_E/£’ 0 )~“~^ log(£ / £ o) 
(Landau et al. 1986; Massaro et al. 2004), where the parameter a 
is the spectral slope at the energy Eg and the parameter ft measures 
the curvature around the peak. We fixed the reference energy Eg to 
300 MeV. The fit with a LP results in a TS = 554, with spectral 
parameters a = 2.35 ± 0.09 and ft = 0.06 ± 0.04. Applying a 
likelihood ratio test to check the PL model (null hypothesis) against 
the LP model (alternative hypothesis), the PL spectral model is 
favoured, indicating that no significant curvature was observed in 
the average y-ray spectrum. 


2 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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Figure 1. Two-month integrated flux ( E > 100 MeV) light curve of PKS 
2123—463 obtained from 2008 August 4 to 2012 January 4 (MJD 54682- 
55930). The last bin covers one month of observation. Arrows refer to 2 a 
upper limits on the source flux. Upper limits are computed when TS < 10. 


Fig. 1 shows the y-ray light curve of the first 41 months of Fermi 
observations of PKS 2123—463 built using two-month time bins, 
with the exception of the final (one-month) data point. For each 
time bin the photon index was frozen to the value resulting from 
the likelihood analysis over the entire period. If TS < 10 2 a upper 
limits were calculated instead. The systematic uncertainty in the 
flux is energy dependent: it amounts to 10 per cent at 100 MeV, 
decreasing to 5 per cent at 560 MeV and increasing to 10 per cent 
above 10 GeV (Ackermann et al. 2012). 

As shown in Fig. 1, PKS 2123—463 was in alow y-ray state (0. 1 — 
100 GeV flux < 5 x 10 8 photons cm 2 s *) for the first 2.5 yr. In 
particular, the source was not detected with two-month time bin 
during the period 2009 August-2010 May. A first increase of the 
y-ray flux was observed in 2011 February-March, and subse- 
quently a flaring activity in 2011 December. We show a light curve 
focused on the period of high activity (2011 December 4-2012 
January 4; MJD 55899-55930) with 2-d or 4-d time bins (Fig. 2). 
The peak of the emission was observed between December 13 
15:43 ut and December 14 15:43 ut, with a flux of (128 ± 23) x 
10~ 8 photons cm~ 2 s _1 in the 0.1-100 GeV energy range, a factor 
of ~25 higher with respect to the average y -ray flux observed during 
2008-20 1 1 . The corresponding observed isotropic y -ray luminosity 
peak in the 0.1-100 GeV energy range is 8.9 ± 0.8 x 10 48 erg s _1 
(assuming z = 1.46 ± 0.05, see Section 6), comparable to the val- 
ues reached by the most powerful FSRQs in flaring activity (e.g. 
Ackermann et al. 2010; Orienti et al. 2012). Leaving the photon 
index free to vary during the period of high activity and Eg fixed to 
670 MeV the fit results in a photon index T = 2.26 ± 0.06, show- 
ing a moderate harder when brighter behaviour already observed in 
other FSRQs (Abdo et al. 2010b). Replacing in the same period the 
PL with a LP, fixing the reference energy Eg to 300 MeV, we obtain 
spectral parameters a = 2.09 ±0.13 and ft = 0.10 ± 0.04. We 
used a likelihood ratio test to check the PL model (null hypothesis) 
against the LP model (alternative hypothesis). These values may be 
compared, following Nolan et al. (2012), by evaluating the curva- 
ture test statistic TS curve = 2[log(like) LP — log(like) PL ] = 3(~1.7cr); 
thus no significant curvature was observed in the y-ray spectrum 
also during the high activity period. 
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Figure 2. Integrated flux (E > 100 MeV) light curve of PKS 2123—463 
obtained from 2011 December 4 to 2012 January 4 (MJD 55899-55930) 
with 2-d or 4-d time bins. Arrows refer to 2a upper limits on the source flux. 
Upper limits are computed when TS < 10. 

Finally, using the gtsrcprob tool, which calculates the proba- 
bilities an event belongs to each of the sources in the model, we 
estimated that the highest energy photon emitted by PKS 2123—463 
was observed with a probability of 87 per cent at distance of 0?35 
from the source with an energy of 92 GeV. For a single high-energy 
photon (> 1 GeV), sky location can only be determined within ~0?8 
(95 per cent confidence radius). 

3 SWIFT DATA: ANALYSIS AND RESULTS 

The Swift satellite (Gehrels et al. 2004) performed two observa- 
tions of PKS 2123—463 in 2011 December triggered by the y-ray 
flare detected by Fermi-LKT. As a comparison we also analysed 
two additional observations performed during 2010-2011. The ob- 
servations were performed with all three on-board instruments: the 
X-ray Telescope (XRT; Burrows et al. 2005, 0.2-10.0 keV), the 
UVOT (Roming et al. 2005, 170-600 nm) and the Burst Alert Tele- 
scope (BAT; Barthelmy et al. 2005, 15-150 keV). 

The hard X-ray flux of this source is below the sensitivity of the 
BAT instrument for the short exposure time of these observations. 
Moreover, the source was not present in the Swift BAT 58-month 
hard X-ray catalogue (Baumgartner et al. 2010) and the 54-month 
Palermo BAT catalogue (Cusumano et al. 2010). 

3.1 Swift /XRT 

The XRT data were processed with standard procedures 
(xrtpipeline vO . 12 . 6), filtering and screening criteria by using 
the heasoft package (v6.11). 3 The data were collected in photon 
counting mode in all observations, and only XRT event grades 
0-12 were selected. The source count rate was low 
(<0.5 counts s~'); thus pile-up correction was not required. Source 
events were extracted from a circular region with a radius of 20 pixel 
(1 pixel ~ 2.36arcsec), while background events were extracted 
from a circular region with a radius of 50 pixel away from the source 
region. Ancillary response files were generated with xrtmkarf , and 

3 http://heasarc.nasa.gov/lheasoft/ 


account for different extraction regions, vignetting and PSF cor- 
rections. We used the spectral redistribution matrices v013 in the 
calibration data base maintained by heasarc. 

Considering the low number of photons collected (<200 counts) 
the spectra were rebinned with a minimum of 1 count per bin and the 
Cash statistic (Cash 1979) was applied. We fitted the spectrum with 
an absorbed PL using the photoelectric absorption model tbabs 
(Wilms, Allen & McCray 2000), with a neutral hydrogen column 
fixed to its Galactic value (2.34 x 10 20 cm~ 2 ; Kalberla et al. 2005). 
The fit results are reported in Table 1 and Fig. 4. 

A marginal increase (1.5cr) of the X-ray flux with respect to the 
previous Swift /XRT observations was observed on 201 1 December 
15, soon after the y-ray flare. The source remained at a similar flux 
level in X-rays on 201 1 December 19. 

3.2 Swi/if/UVOT 

UVOT UV/optical imaging was obtained during all four Swift ob- 
servations of PKS 2123—463. The photometry was carried out on 
pipeline processed sky images downloaded from the Swift data cen- 
tre, 4 following the standard UVOT procedure (Poole et al. 2008). 
Source photometric measurements were extracted from the UVOT 
imaging data using the tool uvotmaghist (vl.l) with a circular 
source extraction region that ranged from 3.5 to 5 arcsec radius to 
maximize the signal-to-noise ratio. In order to remain compatible 
with the effective area calibrations, which are based on 5 arcsec 
aperture photometry (Poole et al. 2008), an aperture correction was 
applied where necessary. This correction was at maximum 5-6 per 
cent of the flux, depending on the filter. The UVOT photometry is 
presented in Table 2 and Fig. 4. Contemporaneous UVOT observa- 
tions in 2011 December 15 found PKS 2123—463 about 0.5 mag 
brighter in v, u and w\ bands, 0.7 mag in ml band, and about 
0.8 mag in wl band compared to the UVOT observation performed 
on 2011 July 3. 

4 GROND DATA 

On 2011 December 18, at 01:21 ut, PKS 2123—463 was observed 
with the GROND (Greiner et al. 2008) mounted on the MPG/ESO 
2.2-m telescope at La Silla, Chile. Preliminary results have been 
reported in Rau et al. (2011). GROND is a seven-channel imager 
that observes in four optical and three near-infrared (near-IR) chan- 
nels simultaneously. The data were reduced and analysed with the 
standard tools and methods described in Rruhler et al. (2008). Cal- 
ibration was performed against a Sloan Digital Sky Survey (SDSS) 
standard star field (g'/i'z!) and against selected Two Micron All 
Sky Survey (2MASS) stars (Skrutskie et al. 2006) (J HKs). This 
resulted in lcr accuracies of 0.05 mag (g'/i'z!, J, H) and 0.07 mag 
(Kg). All magnitudes are corrected for Galactic foreground ex- 
tinction of E(B — V) = 0.030 mag (Schlegel et al. 1998) and are 
summarized in Table 3. 

5 RADIO DATA 

PKS 2 1 23—463 is monitored by the Tracking Active galactic Nuclei 
with Austral Milliarcsecond Interferometry (TANAMI) programme 
(Ojha et al. 2010) at a number of radio frequencies and resolutions 
with the ATCA and the Ceduna facilities. ATCA is observing this 
source every few weeks with ‘snapshot’ observations at frequencies 
from 5.5 through 40 GHz where each frequency is the centre of a 

4 http://www.swift.ac.uk/swift_portal 
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Table 1. Log and fitting results of Swift/X RT observations of PKS 2123—463 using 
a power-law model with Nu fixed to Galactic absorption. 


MJD 

ut date 

Exposure time 
(s) 

r 

Flux (0.3-10 keV)“ 
(xlO” 12 erg cm” 2 s” 1 ) 

55467 

2010-09-28 

2712 

1.47 ±0.30 

1.31 ± 0.46 

55745 

2011-07-03 

2912 

1.61 ± 0.28 

1.20 ±0.32 

55910 

2011-12-15 

4947 

1.68 ±0.18 

1.90 ±0.34 

55914 

2011-12-19 

5521 

1.26 ± 0.18 

1.83 ± 0.27 


“Observed flux. 


Table 2. UVOT photometry of PKS 2123—463. 


ut date 



AB magnitude 0 




uvw2 

uvm2 

uvwl 

u 

b 

V 

2010-09-28 18:29-23:30 

19.78 ± 0.09 

19.34 ±0.10 

19.17 ±0.10 

18.66 ±0.09 

18.61 ±0.14 

18.59 ± 0.20 

2011-07-03 17:15-23:29 

19.83 ± 0.20 

19.28 ±0.14 

18.75 ± 0.11 

18.47 ± 0.12 

18.81 ± 0.25 

18.21 ± 0.26 

2011-12-15 20:45-01:25 

18.97 ± 0.06 

18.58 ± 0.06 

18.24 ±0.05 

17.95 ± 0.06 

- 

17.64 ± 0.06 

2011-12-19 04:20-04:39 

19.17 ± 0.12 

18.80 ±0.12 

18.51 ±0.09 

18.04 ±0.07 

17.77 ± 0.07 

17.82 ±0.15 


“Corrected for Galactic foreground reddening of E(B — V) = 0.03 mag (Schlegel, Finkbeiner & Davis 199S). 


Table 3. GROND photometry of PKS 2123—463. 


ut date 

g' 


AB magnitude 0 
i' z' 

J 

H 

K s 

2011-12-18 01:21 

17.78 ± 0.05 

17.52 ± 0.05 

17.43 ± 0.05 17.24 ± 0.05 

16.75 ± 0.06 

16.53 ±0.06 

16.24 ± 0.09 


“Corrected for Galactic foreground reddening of E(B — V) = 0.03 mag (Schlegel et al. 1998). 


Table 4. ATCA (5.5, 9, 17, 19, 38 and 
40 GHz) and Ceduna (6.7 GHz) observations 
of PKS 2123—463 in 2011. 


ut date 

Frequency (GHz) 

Flux (Jy) 

2011-05-17 

5.5 

0.86 


9.0 

0.77 


17.0 

0.67 


19.0 

0.64 


38.0 

0.45 


40.0 

0.44 

2011-10-14 

5.5 

0.86 


9.0 

0.76 


17.0 

0.60 


19.0 

0.56 


38.0 

0.42 


40.0 

0.41 

2011-11-08 

5.5 

0.85 


9.0 

0.71 


17.0 

0.57 


19.0 

0.53 


38.0 

0.43 


40.0 

0.43 

2011-12-19 

6.7 

0.81 

2012-02-18 

6.7 

0.78 


2 GHz wide band and the fluxes are calibrated against the ATCA 
primary flux calibrator PKS 1934—638 (Stevens et al. 2012). Each 
flux density has a ler uncertainty of 0.01 Jy. The Ceduna radio tele- 
scope in South Australia is monitoring PKS 2123—463 at 6.7 GHz. 
Each flux density has a lcr uncertainty of 0.3 Jy (McCulloch et al. 
2005). No sign of increased activity of the flux density was de- 

©2012 The Authors. MNRAS 427, 893-900 

Monthly Notices of the Royal Astronomical Society © 2012 RAS 


tected between 40 and 5.5 GHz before the y-ray flare (see Fig. 4 
and Table 4). 

In addition, observations of PKS 2123—463 were made after the 
y-ray outburst using the KAT-7 array (the prototype for meerkat) 
in the Karoo. Since the array is still undergoing commissioning 
tests some uncertainties remained about the absolute calibration 
scale. To minimize these uncertainties the observations were taken 
at similar times of days (hence local sidereal time range) in late 201 1 
December to early 2012 January, with 4-5 h durations. The very first 
observation (December 21) failed because the online fringe stop- 
ping was not working properly, but as the array is very small (longest 
baseline 200 m) and the frequency low enough ( 1 .822 GHz) the de- 
lay tracking corrections could be done offline with only marginal 
loss in signal-to-noise ratio. Subsequent observations were done in 
this mode. There were also different problems with various receivers 
over this time so un-coverage and maps were not directly compa- 
rable. The bandpass and gain solutions (on PKS 1934—638) were 
good enough and checks against a nearby phase reference source 
(PKS 2134—470) gave sufficient confidence in the visibilities and 
phases to make simple maps of the inner 0?5 of the primary beam 
and verify that no sources stronger than 100 mjy were present. 

Assuming a flux density for PKS 1934—638 of 13.6 Jy at 
1.822 GHz (as interpolated from ATCA models) the measured flux 
density of PKS 2123—463 is 1.13 Jy on 201 1 December 24, 26, 29, 
and 2012 January 7 and 13. The uncertainties on the absolute flux 
density scale are 5 per cent. As a comparison we extrapolate the flux 
density at 1.82 GHz using the values at 2.7, 0.8 and 0.4 GHz from 
the Parkes Catalogue (1990) and the Molonglo Sky Survey (2008). 
Assuming the spectral index 0.4 derived between 2.7 and 0.4 GHz, 
we obtain a flux density at 1.82 GHz of 1.0 Jy, indicating that no 
significant variation was observed during the y-ray flaring activity 
at this frequency. 
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Figure 3. UV-near-IR SED composed of Svei/f/UVOT observations ob- 
tained on 2011 December 19 and GROND data taken on 201 1 December 
18. The simple power-law model (solid line) and the AGN/galaxy hybrid 
template (dashed line) suggest a photometric redshift of z & 1 -45 (see text 
for details). 

6 NEW PHOTOMETRIC REDSHIFT 
ESTIMATION 

The redshift of PKS 2123—463 has not been established convinc- 
ingly yet. Following the method described in Rau et al. (2012) 
and applied to a sample of 103 Fermi- LAT blazars, we combined 
the GROND photometry with the UVOT photometry from 2011 
December 19 in order to construct a 13-band SED. In order to 
account for source-intrinsic variability between the GROND and 
UVOT observing epochs, the GROND photometry was corrected 
by —0.05 mag in all bands (see Rau et al. 2012, for a detailed 
description of the method). The resulting SED, corrected for the 
Galactic foreground reddening of E(B — V) = 0.03 mag (Schlegel 
et al. 1998), has been fitted with a set of PL models as well as with 
hybrid templates built from normal galaxies and AGN (Salvato et al. 
2009, 2011) using the le phare code (Arnouts et al. 1999; Ilbert et al. 
2006). The best-fitting redshifts for both template libraries are in 
good agreement (see Fig. 3). For the PL model we obtain a 99 per 
cent confidence redshift of z = 1.37^q21 (x 2 = 17, P- = 94.6 s ) and 
for the hybrid models we find z= 1.46 ± 0.05 (y 2 = 23, P- = 99.2). 
This is in rough agreement with the initial spectroscopic redshift of 
Z= 1.67 from Savage & Wright (1981). 

7 SED MODELLING AND CONCLUSIONS 

Beyond the excellent spatial association obtained with the 41 -month 
Fermi-FKt data set, the most secure and distinctive signature for 
firm identification of the y-ray source detected by Fermi- LAT with 
the blazar PKS 2123—463 is the simultaneous increase observed 
in the y-ray, X-ray and optical-UV bands (see Fig. 4; Sections 2 
and 3). 

In order to investigate the physical properties of the source we 
have built a simultaneous SED of the flaring state of PKS 2123—463. 
The Fenni- LAT spectrum was built with data from observations cen- 
tred on 2011 December 10 to 19 (MJD 55905-55916). In addition, 
we included in the SED the GROND and Swift (UVOT and XRT) 
data collected on 201 1 December 18 and 19, respectively. Here, the 
flux suppression, in particular in the UV bands, was corrected as- 


5 The best-fitting templates and redshifts were selected on the basis of a 

simple x 5 2 * method. P z is the integral of the probability distribution function 

f f(z)dz at r p hot ± 0.1(1 + Jphot). which describes the probability that the 

redshift of a source is within 0.1(1 + z) of the best-fitting value. 



T (MJD) 


Figure 4. Light curves of PKS 2123—463 collected in y rays by Fermi- 
LAT (0.1-100GeV flux in units of 10~ 8 photons cm -2 s _1 ), in X-rays 
by Swift/XKT (flux in units of 10 12 ergcnU 2 s _1 ), in UV (filters w 1: 
empty squares, m2: empty triangles and w2: filled circles; flux in units of 
10 12 erg cm 2 s 1 ) and in optical (filters u: empty squares, b: empty trian- 
gles and v: filled circles; flux in units of 10 13 erg cm -2 s _1 ) by Swift/UVOT, 
and radio by ATCA (empty squares: 40 GHz, empty triangles: 17 GHz and 
filled circle: 5.5 GHz) between 2010 June and 2011 December. 

suming a PL spectral shape consistent with the optical-near-IR mea- 
surements (F- a a A 0 - 8 ). The data from ATCA and KAT-7 collected 
on December 19 and 24, respectively, provided information about 
the radio part of the spectrum. The flare centred on ~MJD 55909 
had a variability time-scale of ~2 d, which constrains the size of 
the emitting region during the flare to R' b < 2.1 x 10 16 05 fl /20) cm. 

A ‘blue bump’ accretion disc component is clearly visible in the 
optical/UV data. We modelled these data with a combination of a 
non-thermal synchrotron component and a Shakura-Sunyaev disc 
component (Shakura & Sunyaev 1973). A fit to the disc compo- 
nent allows us to get a rough estimate of the black hole mass of 
M bh « 2 x 10 9 Mq. This is the first black hole mass estimate 
for this source. It is consistent with black hole estimates for other 
high-z FSRQs, obtained in a similar way by Ghisellini et al. (2011). 
Our mass estimate for PKS 2123—463 follows Ghisellini et al. 
(2011) and assumes that the innermost stable circular orbit (ISCO) 
is Rdisc.isco = 6 R g , as one would expect for a Schwarzschild (i.e. 
non-rotating) black hole. If the jet is produced by the Blandford- 
Znajek (Blandford & Znajek 1977) or Blandford-Payne (Blandford 
& Payne 1982) mechanisms this requires a non-zero black hole 
spin, and one expects Raise, isco to be smaller or larger, depend- 
ing on whether the spin is retrograde or prograde (e.g. Garofalo, 
Evans & Sambruna 2010). Due to the uncertainty in spin, this mass 
estimate can be considered to have considerable uncertainty. 

We modelled the portion of the SED from X-rays to y rays 
assuming emission from a relativistic jet with mechanisms of syn- 
chrotron self-Compton (SSC) and Compton scattering of a dust 
torus external to the jet [external Compton (EC)-dust]. The de- 
scription of the model can be found in Finke, Dermer & Bottcher 
(2008) and Dermer et al. (2009). The synchrotron component con- 
sidered is self-absorbed below ~10 u Hz. Correlations of y- ray 
and optical flares with radio light curves and rotations of optical 
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polarization angles in low-synchrotron-peaked blazars seem to indi- 
cate that the y-ray/optical-emitting region is outside the broad-line 
region (BLR), where the dust torus is the likely seed photon source 
(e.g. Marscher et al. 2010). From the disc luminosity obtained by 
the fit of the blue bump, L disc = 1.8 x 10 46 erg s ', we estimate an 
associated BLR radius Rblr = 4.4 x 10 17 cm, based on the relation 
between disc luminosity and Rblr determined from reverberation 
mapping campaigns (e.g. Kaspi et al. 2005; Ghisellini & Tavecchio 
2008). To minimize the scattered BLR contribution, we placed the 
emitting region at r > Rblr- Here the primary seed photon source is 
the dust torus, which was simulated as a one-dimensional ring with 
radius R dllst aligned orthogonal to the jet, emitting as a blackbody 
with temperature r dust and luminosity L dust . 

The model fit to the broad-band SED can be seen in Fig. 5 and the 
parameters can be found in Table 5. A description of the parameters 



Figure 5. Spectral energy distribution data (circles and squares) and model 
fit (solid curve) of PKS 2123—463 with the model components shown as 
dashed curves. The data points were collected by GROND (201 1 December 

18) , Swift (UVOT and XRT; 2011 December 19) and Fermi - LAT (2011 
December 10-19), together with radio data from ATCA (2011 December 

19) and KAT-7 (2011 December 24). 


Table 5. Model parameters for the SED shown in Fig. 5. 


Redshift 

z 

1.46 

Bulk Lorentz factor 

r 

20 

Doppler factor 

<$D 

20 

Magnetic field 

B 

0.8 G 

Variability time-scale 

tv 

1.7 x 10 5 s 

Comoving radius of blob 

K 

4.1 x 10 16 cm 

Jet height 

r 

l.Ox 10 18 cm 

Low-energy electron spectral index 

Pi 

2.0 

High-energy electron spectral index 

P2 

3.8 

Minimum electron Lorentz factor 

Ymin 

4.0 

Break electron Lorentz factor 

Ybrk 

1.0 x 10 3 

Maximum electron Lorentz factor 

Ymax 

1.0 x 10 5 

Disc luminosity 

^disc 

1.8 x 10 46 erg s - 

Black hole mass 

Mbh 

2 x 10 9 M o 

Accretion efficiency 

V 

1/12 

Gravitational radius 


1.76x 10 14 cm 

Inner disc radius 

^discJSCO 

6R g 

Outer disc radius 

^disc_max 

10 *R g 

Dust torus luminosity 

f'dust 

1.5 x 10 45 ergs -1 

Dust torus temperature 

T dust 

1.7 x 10 3 K 

Dust torus radius 

^dust 

3.2 x 10 18 cm 

Jet power in magnetic field 

Pj,B 

3.3 x 10 45 erg s - 

Jet power in electrons 

Pj> par 

1.5 x 10 45 erg s _ 
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can be found in Dermer et al. (2009). The dust parameters were 
chosen so that R d ust is roughly consistent with the sublimation radius 
(Nenkova et al. 2008). The EC-BLR component was calculated 
assuming that the seed photons are from Her and have a luminosity 
of 0. 1 T d ;, L . . 

The electron distribution used, a broken PL with index p\ = 2.0 
below the break at y brk and p 2 = 3.8 for y brk < y', is consistent 
with particles injected with index 2.8, and emission taking place 
in the fast-cooling regime (e.g. Bottcher & Dermer 2002). That is, 
particles are injected between y brk and )/ ax with a cooling electron 
Lorentz factor 

3 m e c 2 

Xcool = ~ ~ 

4ccr T te SC u tat , 

where u tot is the total energy density in the blob frame, which in the 
case of our model fit is dominated by the external energy density. 
In this case y coo i is associated with y mbl , since it is in the fast- 
cooling regime. Also note that in this model fit the magnetic field 
and electrons are nearly in equipartition. Jet powers were calculated 
assuming a two-sided jet. 

The Compton dominance for PKS 2123—463, i.e. the ratio of the 
peak luminosities of the Compton and synchrotron components, is 
«50, which is a rather standard value for powerful blazars (e.g. 
Ghisellini et al. 201 1). A large disc luminosity was estimated from 
the UVOT data, as expected for a powerful FSRQ, withaL d ; sc /L Edd = 
0.2 in agreement with the blazar divide proposed by Ghisellini et al. 
(2009) as a result of the changing of the accretion mode. 

Variability is common in y-ray blazars and provides a power- 
ful tool to associate them definitively with objects known at other 
wavelengths and to study the emission mechanisms at work. The 
combination of deep and fairly uniform exposure over ~3 h, very 
good angular resolution, and stable response of the Fermi-LAT is 
producing the most sensitive, best-resolved survey of the y-ray sky. 
On the other hand, those cases where there is a decrease in the ac- 
tivity, and thus of the significance of detection, can lead to a more 
complex identification process for a y-ray source. When combined 
with simultaneous ground- and space-based multifrequency obser- 
vations, the Fermi-LAT achieves its full capability for the identifi- 
cation of the y-ray sources with counterparts at lower energies and 
the knowledge of their emission processes, as reported here for the 
high-z, Compton dominated FSRQ PKS 2123—463. 
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